Computational Biology and Chemistry

Abstract Background: The stem cell factor receptor or c-Kit is a type III receptor tyrosine kinase, activated by its ligand Stem cell factor (SCF). Up on activation, c-kit induces signaling pathways that regulates blood cell proliferation, survival, differentiation, and migration. Several studies reported that c-Kit/SCF signaling, contributes to the development and progression of acute myeloid leukemia (AML) in patients. However, the downstream proteins regulated by c-kit activation and their clinical significance in AML remain poorly explored. Methods: Human Acute megakaryoblastic leukemia (Mo7e) cells, were-stimulated with SCF and global protein expression were profiled using two-dimensional gel electrophoresis coupled with MALDI-TOF and LC-MS/MS. Differentially expressed proteins were functionally characterized and validated using patient data from the TCGA-LAML and matched normal data from GTEx, GEO datasets, and quantitative RT-PCR. Their diagnostic and prognostic significance was assessed using ROC, Cox regression, LASSO, Kaplan Meier survival analyses, and a prognostic nomogram model. Results: Proteomic profiling identified 14 differentially expressed proteins in SCF-stimulated Mo7e cells, which are predicted to involved in cytoskeletal organization, protein folding, metabolism, vesicular trafficking, and translational regulation. Transcriptomic analysis of the TCGA-LAML cohort revealed significant dysregulation of CFL1, CCT8, HSP90B1, MDH2, EIF5A, GSN, and TPI1. Integrated ROC, Cox regression, and LASSO analyses identified CFL1, CCT8, and GSN as the most robust prognostic biomarkers associated with poor overall survival in LAML patients. Their expression patterns were validated in independent GEO datasets and by qRT-PCR in SCF stimulated Mo7e cells. Finally, a three-gene nomogram model was developed and validated to predict the overall survival probability of AML patients at 1-, 3-, and 5-year time points. Conclusions: This study identifies CFL1, CCT8, and GSN as key downstream effectors of c-Kit signaling as prognostic biomarkers for AML. These findings provide mechanistic insights into c-Kit-driven leukemogenesis and establish a clinically relevant three-gene signature for AML risk stratification and potential therapeutic targeting.

Aptamers are structured nucleic acid ligands capable of high affinity, high specificity molecular recognition generated using variations of the SELEX (Systematic Evolution of Ligands by Exponential Enrichment) process. However, SELEX often produces sequences that enrich yet may lack binding efficacy. We propose a measure called the Ruggedness Composite Index (RCI) along with a method for computing it, that can be used to distinguish binding-competent ('active') aptamers from weak or non-binding ('inactive') aptamers. Given a set of aptamers, RCI incorporates information on their fragmentation (landscape partitioning), basin entropy (metastable state distribution), cumulative density irregularity (non-uniform occupancy), and structural energy correlation length (structure-energy coupling scale). We test whether secondary-structure folding energy landscape topology distinguishes active from inactive aptamers using a multiscale level set framework across six datasets. Active aptamers show lower RCI values and occupy smoother, funnel-like conformational spaces, while inactive aptamers show higher RCI values, reflecting fragmented, high-entropy landscapes. By contrast, classical thermodynamic features, such as minimum free energy, show limited discrimination between active and inactive aptamers. In all datasets, sequences that exhibit enrichment which is not monotonic but lack specificity exhibit elevated ruggedness, indicating landscape topology can predict non-specific enrichment. These results indicate that folding landscape organization can be used as a predictor of aptamer activity and establish RCI as a simple, mechanistically interpretable measure for improving candidate prioritization, especially in therapeutic aptamer discovery.

This study emphasizes the potential of computational techniques in cancer risk assessment, lighting opportunities for specific and data-driven healthcare solutions. This study examines the use of artificial intelligence (AI), machine learning (ML), and deep learning (DL) approaches to improve cancer risk assessment using a Kaggle dataset. The study uses Java-based ML software to create and evaluate multiple predictive models, taking advantage of its powerful libraries and frameworks for processing and analyzing cancer risk indicators. This work analyzes model performance using 10-fold cross-validation, resulting in reliable generalization and accuracy estimates. Several classification techniques, such as Random Forest (RF) logistic regression (LR), decision trees (DT), Naive Bayes (NB), and Multi-layer perceptron (MLP), are used to assess their efficacy in predicting risk levels for various cancer types. To measure classification effectiveness, key performance metrics such as accuracy, precision, recall, and F1 score are produced, in addition to multi-class confusion matrices. The results show that the RF model is the best classifier for classification, with accuracy of 99.85%, F-measure of 99.80%, precision of 99.80%, and sensitivity of 99.90%. These findings demonstrate the model's ability to effectively estimate cancer risk levels among individuals. of cancer risk estimations, allowing for earlier discovery and more effective medical care.

Melittin exhibits antitumor activity in cervical cancer models, yet the long non-coding RNA (lncRNA) response and associated regulatory networks remain poorly understood. Here, strand-specific RNA-seq data from melittin-treated and untreated U14 murine cervical cancer cells were analyzed to characterize melittin-responsive lncRNAs and explore their potential functional associations. A total of 28,162 lncRNAs were identified, including 27,307 known and 855 novel transcripts. Differential expression analysis revealed 404 differentially expressed lncRNAs (DElncRNAs), comprising 191 upregulated and 213 downregulated lncRNAs, w most of which were predicted to localize to the cytoplasm or nucleus. Cis-target analysis identified 52 neighboring mRNAs as putative targets of 46 DElncRNAs. Functional enrichment highlighted mitochondrial electron transfer and redox-related processes, including the mitochondrial electron transfer flavoprotein complex, electron-transferring-flavoprotein dehydrogenase activity, ubiquinone binding, and quinone binding. In parallel, melittin induced mitochondrial membrane depolarization and increased intracellular reactive oxygen species accumulation in U14 cells. Co-expression analysis further identified 138 lncRNAs co-expressed with 161 mRNAs, which were enriched in chromatin remodeling, DNA replication, and DNA repair. EdU incorporation decreased with increasing melittin concentrations, indicating suppression of DNA synthesis and proliferative activity. RT-qPCR analysis confirmed the expression trends of selected DElncRNAs. Collectively, these findings demonstrate extensive remodeling of the lncRNA landscape in melittin-treated U14 cells and suggest that melittin-responsive lncRNA-mRNA networks are associated with mitochondrial redox disruption and impaired DNA synthesis. This study provides a transcriptomic framework for identifying candidate lncRNA-mRNA regulatory axes underlying the antitumor response to melittin.

A protein whose removal or deletion causes significant disruption or collapse of a protein-protein interaction (PPI) network is referred to as a vulnerable protein. Such proteins may serve as valuable therapeutic or diagnostic targets in disease-associated networks. In this study, two PPI networks were constructed, one for HPV-positive and the other for HPV-negative head and neck squamous cell carcinoma (HNSCC), and the vulnerable proteins of these networks were identified by the node deletion approach. After analyzing the networks, 27 unique vulnerable proteins in HPV-positive and 72 unique vulnerable proteins in HPV-negative HNSCC were identified. Among them, one HPV-positive and seven HPV-negative HNSCC vulnerable proteins were further chosen by integrating multi-omics data. To exploit the vulnerabilities of these proteins, candidate synthetic lethal (SL) partners were predicted whose inhibition may selectively impair tumor survival. Subsequently, drug-gene interaction analysis was performed to identify inhibitors targeting the SL partners of these vulnerable proteins. Notably, in HPV-positive HNSCC, TOP2A, CHEK1, and CHEK2 genes were identified as SL partners of TTN, and their inhibitors were already clinically approved. While in HPV-negative HNSCC, ADA and MMP19 were identified as an SL partner of LMO7; TMEM45B, CDH3, and ELF3 genes were identified as an SL partner of CGN; and ZNF433 was identified as an SL partner of FLNC. However, MMP19, ZNF433, and TMEM45B inhibitors were not reported. Thus, these vulnerable proteins, including their SL partners, provide novel avenues to explore and develop more efficient and precise therapeutic and diagnostic strategies.

Objective White smoke inhalation injury (WSI) causes severe acute lung damage with no specific therapy currently available. Sphingolipid metabolism is implicated in pulmonary inflammation, but its transcriptional regulatory landscape in WSI remains unexplored. This study aimed to identify key sphingolipid metabolism related genes and evaluate their regulatory roles and therapeutic potential in WSI. Methods We established a rat model of WSI and performed integrated bulk RNA sequencing, weighted gene coexpression network analysis (WGCNA), and single-cell RNA sequencing (scRNAseq) to screen for differentially expressed sphingolipid metabolism-related genes (DESRGs). Protein-protein interaction (PPI) network with four centrality algorithms was used to prioritize hub genes. In silico gene knockout and molecular docking were conducted to assess regulatory functions and identify potential drug candidates. Results We identified 22 DESRGs that were predominantly enriched in DNA replication and cell cycle pathways rather than canonical sphingolipid metabolic processes. PPI consensus prioritized three hub genes--Top2a, Ttk, and Ccna2--with Top2a exhibiting the highest expression in epithelial cells and significant downregulation after smoke exposure. ScRNAseq revealed immune cell infiltration and epithelial differentiation trajectories. Virtual knockout showed that Top2a depletion affected the largest transcriptomic fraction (~0.4%) and was enriched in lysosome biogenesis, innate immunity, phagocytosis, and lipid catabolism. Molecular docking identified thalidomide as a high affinity ligand for Top2a (Vina score: -8.5 kcal/mol). Conclusion Our multiomics integrative framework identifies Top2a as a central regulatory hub linking sphingolipid associated inflammation to epithelial responses in WSI, and nominates thalidomide as a potential drug repurposing candidate. These findings provide prioritized targets for future translational investigation.

Archaea are promising chassis organisms in biotechnology, and the accurate annotation of their chromosomal replication origins (oriCs) is the key to unlocking their full potential. However, the existing Ori-Finder 2 web server suffers from low accuracy, slow speed, and limited scalability. In this study, we present Ori-Finder-Arch, an updated web server for high-performance oriC prediction in archaea. This pipeline integrates HMMER-based replication initiation protein (RIP) annotation, refined consensus motif recognition, and GC profile-based DNA unwinding element (DUE) detection. On a benchmark set of experimentally validated oriCs, Ori-Finder-Arch achieved a recall of 95.6% and a precision of 86.0%, substantially outperforming Ori-Finder 2 (62.2% and 63.6%, respectively), while running 4.75 times faster and supporting diverse assembly levels. When applied to the available archaeal assemblies, it successfully annotated 17,472 oriCs. Meanwhile, the web server provides interactive visualizations at different levels. In conclusion, Ori-Finder-Arch offers an efficient, accurate, and user-friendly platform for advanced studies of archaeal DNA replication initiation and synthetic biology applications, and is freely available at https://tubic.org/Ori-Finder-Arch/ and https://tubic.tju.edu.cn/Ori-Finder-Arch/.

Background: RNA molecules form thermodynamic ensembles, but interpretation often requires a single representative structure. Existing base-pair centroid estimators assess agreement at the level of individual base pairs and do not directly target nesting depth along the sequence. Methods: We introduce Mountain Centroid, which minimizes expected squared mountain-profile distance, and derive dynamic programming algorithms with and without RNA pairing constraints. We also combine the Mountain Centroid objective with the base-pair centroid gain. Results: Across 21,254 RNAStrAlign sequences, Mountain Centroid had lower median normalized mean squared mountain distance (NMSMD) than minimum-free-energy (MFE) and base-pair centroid ({gamma} = 1) structures, whereas its median base-pair F1 was lower. Imposing RNA pairing constraints improved base-pair F1 for 59.35% of sequences and reduced it for 3.58%. At an illustrative weight, the combined objective had median base-pair F1 similar to MFE while retaining lower median NMSMD than MFE and all tested {gamma}-centroid settings. Conclusions: Mountain Centroid represents an RNA structural ensemble with a single secondary structure that reflects how nesting depth varies across nucleotide positions. Combining mountain-profile and individual-base-pair criteria allows their relative contributions to be varied.

Recent advances in artificial intelligence have accelerated the discovery of bioactive peptides by enabling computational exploration of the vast peptide sequence space. However, existing peptide generation approaches generally rely on either distribution-learning models, which generate biologically realistic sequences but do not consistently optimize functional activity, or optimization-based methods, which maximize prediction confidence while often deviating from the underlying distribution of experimentally validated peptides. To address this limitation, a two-phase generative-evolutionary framework is proposed that integrates distribution learning with evolutionary optimization. In the first phase, Variational Autoencoders (VAE), Autoregressive Transformers (ART), and Token Diffusion Transformers (TDT) are used to generate biologically plausible seed peptides. In the second phase, these peptides were used as initial seed for Hill Climbing optimization procedure that iteratively improves fitness function score. The proposed two-phase framework was evaluated using a dataset of experimentally validated IL-2-inducing peptides. Evaluation using independent IL-2 prediction models showed that Autoregressive Transformer combined with Hill Climbing achieved the best overall performance, achieving the mean IL-2 induction confidence score of 0.96 while reducing KL divergence from 2.26 for standalone Hill Climbing to 0.75. A case study on an independent IL-13 inducing peptide dataset showed similar trends, with ART initialized Hill Climbing achieving the mean IL-13 induction score of 0.99 while reducing KL divergence from 1.76 to 0.59. Overall, the framework provides a generalizable approach for balancing functional optimization and distributional realism and can be applied to peptide discovery and data augmentation in imbalanced biological datasets thereby generating high confidence peptides for wet lab validation. HighlightsO_LIProposed a two-phase framework for bioactive peptide generation with potential to address class imbalance in peptide classification tasks. C_LIO_LIPerformed a systematic comparison of distribution-learning and optimization-based approaches for peptide generation. C_LIO_LICombined distribution-learning models for sequence generation with optimization algorithms for improving peptide functional properties. C_LIO_LIDemonstrated the applicability of the proposed framework across multiple bioactive peptide datasets. C_LI

The attention mechanism offers the possibility for data-driven discovery of biological principles. However, for important protein families such as human olfactory receptors, the extent to which attention can associate with biologically meaningful key regions lacks systematic validation. In this study, using human olfactory receptors (ORs) as a model, we constructed CrossVOI, a VOC-OR interaction prediction framework based on protein language models and cross-attention, achieving predictive performance superior to existing methods. Furthermore, we systematically analyzed the attention distributions of CrossVOI and found that attention not only focused on ligand-binding interfaces and evolutionarily conserved sites, but also to some extent identified certain dynamically regulated regions. In summary, we propose CrossVOI, currently the best-performing framework for VOC-OR interaction prediction, and analyze the interpretability of the attention mechanism for human ORs. This study provides insights into the interpretability of protein function prediction methods and is expected to contribute to the exploration of attention mechanisms in biological mechanisms, and provide assistance for large-scale screening and mechanistic analysis of olfactory receptors.

Background: Spatial organization is increasingly recognized as a key determinant of tumor-immune interactions in head and neck squamous cell carcinoma (HNSCC). The GSE300147 Xenium spatial transcriptomic resource generated by McCord and colleagues established a framework for mapping spatially coordinated T-cell states in HNSCC. However, how tumor-enriched epithelial immune states relate to metabolic, redox, and stress-adaptive transcript programs remains incompletely defined. Methods: A secondary, data-driven reanalysis of GSE300147 was performed, focusing on 17 confirmed HNSCC Xenium sections after exclusion of a non-HNSCC ameloblastoma specimen. A total of 1,148,244 cells were analyzed, including 558,867 EpCAM+ tumor-enriched epithelial cells. Tumor-enriched epithelial cells were classified into Hot, Intermediate, and Cold states using a Composite Hotness framework integrating T-cell inflammatory signature score, checkpoint-associated signaling, CD274 expression, IFN/antigen-presentation signature score (IFN/AP), and tumor-immune proximity. Six metabolic ecosystem states, neighborhood profiling, spatial permutation testing, and an integrated Immune-Metabolic-Redox Ecosystem Score (IMRES) were then applied. Results: Immune activation was spatially heterogeneous across HNSCC sections. Immune-hot tumor-enriched epithelial regions showed not only inflammatory, checkpoint-associated, and antigen-presentation signature scores, but also coordinated metabolic, oxidative-redox, and stress-response transcript programs. IMRES, derived from available immune, metabolic, redox, and stress-response transcript components represented in the Xenium panel, increased progressively from Cold to Intermediate to Hot tumor-enriched epithelial states and was associated with NFE2L2, GDF15, HLA-DRA, CD274, KEAP1, and MDM2. Integrating IMRES with Composite Hotness identified a distinct Hot+IMREShigh ecosystem comprising 106,874 tumor-enriched epithelial cells. This state showed the strongest immune-active and stress-adaptive features and was positioned closer to immune populations than expected by random assignment. An alternative rank-based robustness analysis reproduced the IMRES-associated ecosystem axis and correlated with the original module-based score (Spearman r = 0.597). Conclusions: This secondary reanalysis extends the original spatial T-cell framework by defining a complementary tumor-centered immune-metabolic-redox ecosystem in HNSCC. IMRES provides a transcript-derived framework for identifying Hot+IMREShigh neighborhoods where immune activation, checkpoint signaling, metabolic remodeling, and stress adaptation converge, providing a hypothesis-generating framework for studying immune resistance and therapeutic vulnerability.

1. Plants produce a great chemodiversity, which is the diversity of specialized metabolites (SMs). These SMs are produced in complex metabolic pathways and play an important role in inter-species interactions. There are numerous hypotheses about the evolutionary processes which brought about and maintain chemodiversity. Some have been partially tested in lab and field studies. However, some of their assumptions and predictions are better tested by quantitative modeling, and so far no quantitative model has investigated the role of metabolic pathways. 2. To close this gap, we developed an individual-based model for metabolic pathway evolution. It models enzymes creating metabolites with various modifications. Enzymes undergo inheritance and mutation. We used the model to compare the screening and interaction diversity hypotheses. 3. The screening hypothesis predicts promiscuous enzymes, genetic drift, the presence of many non-beneficial metabolites, and high metabolite richness. The interaction diversity hypothesis predicts specialized enzymes, selection, the almost exclusive presence of beneficial metabolites, and situation- dependent metabolite richness. We found that the patterns predicted by the screening hypothesis did not occur, while those predicted by the interaction diversity hypothesis did. 4. This provides reason to favor the interaction diversity hypothesis over the screening hypothesis when connecting empirical results to their evolutionary context

Immune checkpoint blockade has revolutionized cancer therapy, but the therapeutic efficacy is limited. Clinical trials on blockade of newly identified immune checkpoints didn't show promising result, suggesting that it might be insufficient to understand the function of immune checkpoints in cancer merely in the context of immunity. Here, we found mutually exclusive expression patterns of the immune checkpoint VISTA (or VSIR) and the neural stemness factor SETDB1, an oncoprotein that promotes immunoevasion, in xenograft tumors, suggesting that cells with high VISTA expression represents a differentiated, and hence, less or non-malignant state in tumor. Non-neural differentiation factors HHEX, MYOD1 and PPARG promote, whereas oncoproteins KRAS (and the mutant KRAS(G12D)) and SOX2, both being embryonic neural factors, repress VISTA expression. This tendency can be inferred from the finding that neural stemness is the core property of cancer cell. Manipulated expression of VISTA in cancer cells generated no significant effect on cell tumorigenicity and differentiation state, but led to change in cell morphology and actin cytoskeleton. Mechanistically, VISTA regulates a key cytoskeleton regulator, WASF2, leading to the change in cell morphology, which might interfere with signal transduction of immune response. The results suggest that 1) high expression of a protein in tumor might represent a less or non-malignant state, targeting of which would leave malignant cells intact, and consequently, leading to weak or even no therapeutic efficacy, a key factor worth considering for target selection; 2) immune checkpoints might play other roles in cells that interfere with regulation of anti-tumor immunity.

Lung cancer continues to be the leading cause of cancer-related mortality globally, with non-small cell lung cancer (NSCLC) representing the most prevalent subtype. Tumor hypoxia is a characteristic feature of the neoplastic microenvironment in NSCLC, facilitating tumor progression and conferring resistance to oxidative stress through the stabilization of hypoxia-inducible factor-1 alpha (HIF-1). In this study, we investigated the combined anticancer effects of baicalein (Bai), a natural flavonoid, and thermal-cycling stimulation (TCS), a physical treatment that minimizes damage to normal cells, under cobalt (II) chloride (CoCl2)-induced hypoxic conditions in NSCLC. In A549 NSCLC cells, the combination of Bai and TCS significantly decreased cell viability and induced apoptosis, while exhibiting minimal cytotoxicity on IMR-90 normal human lung fibroblast cells. On a mechanistic level, this combined treatment suppressed the expression of HIF-1 and superoxide dismutase 2 (SOD2) proteins, elevated intracellular reactive oxygen species (ROS) levels, and impaired DNA repair capability by downregulating MutT homolog 1 (MTH1) protein expression. Additionally, disruption of mitochondrial membrane potential and increased poly (ADP-ribose) polymerase (PARP) cleavage further confirmed the induction of apoptosis. These findings indicate that combining Bai with TCS offers a promising synergistic approach to treating NSCLC under hypoxic conditions.

Chemoradiotherapy (CRT) is the standard-of-care therapy for many solid malignancies, yet predictive biomarkers of treatment response remain limited. We identified a germline single nucleotide polymorphism (SNP) in an intrinsically disordered region of the lysine demethylase KDM3C/JMJD1C (p.S464T) that is associated with CRT outcomes in locally advanced rectal cancers (LARC) and head and neck squamous cell carcinoma (LA-HNSCC). In silico modeling with AlphaFold predicted S464T substitution influenced interaction between phosphorylated KDM3C and RNF8 FHA domain. In cellular models, conversion of S464 to T464 increased sensitivity to DNA-damaging agents. S464T substitution impaired damage-induced MDC1-RAP80 signaling and downstream RAP80-BRCA1 colocalization. SNP carrying cells impaired DNA repair causing genotoxic stress that is associated with increased cGAS-cGAMP innate immune signaling and increased apoptosis. Population analyses with the SNP highlighted an increase incidence of UV-induced skin and other cancers, linking inherited variation in the chromatin regulatory gene KDM3C to genome instability, cancer risk, and therapeutic vulnerability.

Click-qPCR is a browser-based application for relative qPCR analysis that requires a tidy-format CSV file containing four columns: sample, group, gene, and Cq. Preparing this input from qPCR instrument output typically requires manual reformatting and calculation of mean Cq values for technical replicates. To simplify this process, we developed Click-Prep (https://kubo-azu.shinyapps.io/Click-Prep/), an interactive web-based application designed specifically to create Click-qPCR input files. Click-Prep imports CSV, TXT, TSV, and XLS/XLSX files and supports skipping of instrument-generated metadata rows, interactive column mapping, and manual assignment of experimental groups. Users can review technical-replicate measurements, exclude selected rows according to predefined quality-control criteria, and calculate mean Cq values for each sample-group-target combination. Missing or nonnumeric Cq values are flagged for review and must be resolved before the mean is calculated. Click-Prep can also combine compatible formatted CSV files, such as datasets obtained from separate qPCR plates. The resulting dataset is exported as a standardized CSV file containing the four fields required by Click-qPCR. By integrating these operations into a guided browser-based workflow, Click-Prep enables users to prepare Click-qPCR input files rapidly and consistently without programming.

Glioblastoma (GBM) is a highly aggressive brain tumor with an extremely poor 5-year survival rate of 6.9%, largely attributable to the lack of reliable biomarkers. While competing endogenous RNA (ceRNA) and copy number variation (CNV) analyses offer unique biomarker identification potential, current approaches neglect the integration of multiple regulatory mechanisms for biomarker detection. To address this limitation, we applied relational graph convolutional networks (RGCNs) to ceRNA and CNV knowledge graphs through a novel late fusion ensemble architecture. The proposed architecture outperformed baseline models and identified five novel biomarkers, including hsa-miR-196a and hsa-miR-224. Kaplan-Meier survival analysis and Cox regression indicated that the identified genes hold significant prognostic and diagnostic power. The early stratification of the Kaplan-Meier curves indicates the potential these genes hold for patient survival prediction. The results illustrate that a late fusion RGCN ensemble effectively captures complex gene interactions, overcoming limitations of existing models and providing a framework for biomarker discovery. The novel biomarkers serve as prospective targets for future GBM therapeutic development and candidates for non-invasive diagnostic assays.

Background Zinc finger protein 36-like 1 (ZFP36L1) is an AU-rich element-binding RNA-binding protein that regulates post-transcriptional gene expression and has been implicated in tumor progression, cell-cycle regulation, and DNA damage responses. However, its functional role in triple-negative breast cancer (TNBC) remains poorly understood. This study investigated the effects of CRISPR/Cas9-mediated ZFP36L1 knockout on cell proliferation, doxorubicin (DOX) sensitivity, cell-cycle progression, and DNA damage responses in MDA-MB-231 TNBC cells. Methods Wild-type (WT) and CRISPR/Cas9-generated ZFP36L1 knockout (KO) MDA-MB-231 cells were cultured under standard conditions. Cellular proliferation was evaluated by cell counting over three weeks. Cell viability following DOX treatment was determined using the MTT assay, and half-maximal inhibitory concentration (IC50) values were calculated. Cell-cycle distribution was assessed by propidium iodide flow cytometry after 24 h of DOX exposure, while DNA damage was quantified by {gamma}-H2AX flow cytometric analysis. Statistical significance was determined using Student's t-test with P < 0.05 considered significant. Results ZFP36L1 knockout reduced the proliferative capacity of MDA-MB-231 cells compared with WT cells. Both cell lines exhibited dose-dependent decreases in viability following DOX treatment. KO cells demonstrated a higher mean IC50 than WT cells (9.64 vs. 8.40 M), indicating a trend toward reduced DOX sensitivity; however, this difference was not statistically significant (P = 0.569). Flow cytometric analysis revealed enhanced accumulation of KO cells in the S and G2/M phases following DOX treatment, suggesting altered cell-cycle checkpoint regulation. Furthermore, KO cells exhibited elevated basal {gamma}-H2AX expression and greater DOX-induced {gamma}-H2AX accumulation than WT cells, indicating increased DNA damage and impaired maintenance of genomic stability. Conclusions CRISPR/Cas9-mediated loss of ZFP36L1 suppresses proliferation, alters cell-cycle checkpoint dynamics, and enhances DNA damage accumulation in MDA-MB-231 TNBC cells. These findings indicate that ZFP36L1 plays a context-dependent role in regulating genomic stability and cellular responses to genotoxic stress, highlighting its potential as a biomarker and therapeutic target in triple-negative breast cancer.

Although most human protein coding genes have functional annotations in databases, such as GeneCards, many remain poorly characterized. To address this gap, computational tools can be leveraged to predict the functional roles of under-annotated genes by extracting patterns from complex biological networks. Here we introduce Brain-for-Biotech (BfBio), a framework designed to identify genes important for vascular endothelial cells (EC), which are crucial cells for vessel formation (angiogenesis), vascular homeostasis, hemostasis and blood/tissue barrier function but also critical mediators of immunity and cancer progression. BfBio utilizes a Personalized PageRank (PPR) algorithm on an integrated network of different omics datasets and publicly available gene-gene/protein-protein interaction databases. In this study, we apply the predictive capabilities of BfBio to infer angiogenic stalk cell phenotype function in genes for which this function was not known before. By leveraging a set of genes characterizing the stalk cell cluster in lung tumor EC models previously identified, we have achieved a high Area Under Receiver Operative Characteristic (AUC-ROC) performance (0.837). Enrichment analysis, coupled with a text mining application, further confirmed that among the 49 predicted genes four of them were poorly characterized yet possessed biologically relevant properties and were linked to cancer, thereby validating BfBio as a robust tool for prioritizing novel therapeutic targets in vascular biology.

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.